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ABSTRACT 

We investigate the ability of high resolution X-ray spectroscopy to directly probe the grain composition of the 
interstellar medium. Using iron K-edge experimental data of likely ISM dust candidates taken at the National 
Synchrotron Light Source at Brookhaven National Laboratory and the Advanced Photon Source at Argonne Na- 
tional Laboratory, we explore the prospects for determining the chemical composition of astrophysical dust and 
discuss a technique for doing so. Focusing on the capabilities of the Astro E2 XRS micro-calorimeters, we assess 
the limiting effects of spectral resolution and noise for detecting significant X-ray absorption fine structure signal 
in astrophysical environments containing dust. We find that given sufficient signal, the resolution of the XRS will 
allow us to definitively distinguish gas from dust phase absorption, and certain chemical compositions. 

Subject headings: dust: XAFS, ISM, composition - ISM: dust composition - technique: solid state - X-rays: 
ISM, dust 
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1. INTRODUCTION 

UV and IR studies have shown that the heavier elements such 
as C, N, O, Mg, Si, and Fe commonly condense to form solid 
particles (Draine 2003; Li & Greenberg 2003; Whittet 2003, 
and references therein for reviews). Combined with planetary 
studies of meteorites and dust, common forms of interstellar 
grains include ice (H>0), graphite (C), silicates (SiOa; ferrosil- 
icate, FeSiOs; fayalite, Fe2Si04; enstatite, MgSiO^; and olivine 
Mg2Si04), and various iron species (pure Fe; wiistite, FeO; 
hematite, FeiOs; and magnetite, FesO/t). UV observational ev- 
idence that a fraction of the heavier elements can be found in 
grains is inferred from the reduction of certain elements from 
the gas phase, or the depletion factor, defined to be the ratio (in 
gas-phase abundances) of the amount expected-to-measured. 
The presence of grains is also inferred from studies in the ra- 
dio and optical. IR emission spectra, while able to distinguish 
some complex molecules through probes of vibrational modes 
(through excitation of phonons rather than electrons, as e.g. in 
the X-rays) are limited by the grain types that it can sensitively 
probe (i.e. mostly polycyclic aromatic hydrocarbon or PAHS, 
graphites, and certain silicates). 

Dust particles smaller than 10 /mi are partially transparent 
to X-rays so that the measured absorption in this energy band 
should be sensitive to all atoms in both gas and solid phase. 
High resolution X-ray spectroscopy therefore provides a pow- 
erful and unique tool for determining the state and composi- 
tion of the material in the interstellar medium (ISM). For inter- 
stellar grains, we would detect oscillatory modulations near the 
photoelectric absorption edge known as X-ray absorption fine 
structure (XAFS). The XAFS is a condensed state modulation 
of the atomic absorption cross-section, thus can distinguish ab- 
sorption by the gaseous and the condensed states. When a pho- 
ton excites a deep-core electron into a high-lying, unoccupied 
electronic state, the wave function of the outwardly propagating 
photoelectron interferes with the portion that scatters from the 



surrounding atoms. This interference produces an oscillatory 
fine structure which is characteristic of the chemical species of 
the absorber. 

The first X-ray study identifying absorption features with the 
ISM were by Schattenburg & Canizares (1986) with the Ein- 
stein Focal Plane Crystal Spectrometer pointings of the Crab 
Nebula. To date, several ISM studies using highly absorbed 
bright X-ray binaries with the Chandra HETGS and LETGS in 
the soft (< 1 keV) X-ray have identified absorption features as- 
sociated with ISM material consisting of atomic and/or molec- 
ular forms of oxygen, neon and iron (e.g. Paerels et al. 2001; 
Takei et al. 2002; Schulz et al. 2002; Juett et al. 2004), but 
with no definitive indication that the absorption features had 
to come from grains of specific chemical makeup. Lee et al. 
(2001, 2002) have claimed direct detections of iron and silicate 
grains in Chandra studies of several black hole systems in the 
soft X-rays. In particular, strong modulations in the form of 
XAFS seen near the Si K edge at ~ 1.84keV (-6.7A) led Lee 
et al. (2002) to conclude that they were associated with sili- 
cate grains in the environment of the source or along the line- 
of-sight. Unfortunately, this was reported only as an exciting 
possibility of a direct detection of XAFS, as a consequence of 
the lack of statistics and because there are Si02 and polysili- 
con in the ACIS detectors. Recently, XAFS detection at Si K 
as well as the K-edges of sulfur and magnesium, attributed to 
ISM dust have been reported by Ueda et al. (2004) in the spec- 
tra of the highly absorbed X-ray binaries GX 13+1, GX 5-1, 
and GX 340+0, although Schulz (private communication) cau- 
tion that some of these features may be instrumental in origin. 
Chandra X-ray spectra have also revealed detailed edge struc- 
ture at > 0.7 keV near the Fe L photoelectric edge which have 
been attributed to dust (e.g. Paerels et al. 2001; inferred in Lee 
et al. 2001). 

While progress is continuing in ISM studies using the Chan- 
dra gratings in the soft X-ray, we consider study of XAFS in 
the hard X-ray to be a more promising mechanism for reveal- 
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ing the chemical species of the condensed state material in the 
ISM. Because our better prospects for X-ray grain studies cur- 
rently rely on pointings towards bright X-ray binaries (XRBs) 
with large amounts of absorption, soft X-ray ISM studies of- 
ten have the added worry of needing to separate any solid state 
features from the myriad of ionized absorption and/or emission 
features originating from the XRB environment (e.g. the ac- 
cretion disk corona or the surrounding ionized plasma). Grain 
studies in the ~7.1 keV (1.75 A) Fe K regime therefore allow 
a better search for XAFS, due to the reduced number of spec- 
tral features from atomic processes in highly ionized plasma 
occurring in this spectral region. Furthermore, the chemical 
shift associated with the oxide species is more pronounced for 
the K-edge spectra than for L-edge spectra, thus more likely 
to be resolved using extant instrumentation. In this article, we 
present K-edge XAFS data of various iron compounds taken at 
the National Synchrotron Light Source (NSLS) at Brookhaven 
National Laboratory (BNL) and at the Advanced Photon Source 
(APS) at Argonne National Laboratory (ANL). Using tech- 
niques developed for the analysis of synchrotron data, we dis- 
cuss the prospect of studying the composition of iron grains in 
the ISM with the XRS calorimeter on Astro E2. Soft X-ray 
XAFS synchrotron studies, and its relevance to studies using 
Chandra, XMM-Newton and future X-ray missions will be pre- 
sented in a later paper - These measurements, together with the 
FeK measurements presented in this paper will eventually be 
jointly analyzed to better constrain the true crystalline structure 
of the material. Ultimately, combining laboratory and space- 
based studies of these features as well as combining hard and 
soft X-ray measurements to compare all available absorption 
edges in the ISM candidate materials will allow us to investi- 
gate the content of the ISM by distinguishing between atomic 
and condensed state absorption in the spectra from astrophysi- 
cal sources. 

2. THE THEORY BEHIND MEASURING XAFS 

We briefly present theory relevant for the interpretation, pro- 
cessing, and analysis of XAFS data. For detailed discussions of 
XAFS theory and practice, see e.g. Rehr & Albers (2000), Kon- 
ingsberger & Prins (1988), Stohr (1996), and/or Kruegel (2003) 
and references therein. See also Woo (1995), Woo et al. (1997), 
and Forrey et al. (1998) for early discussions of XAFS in the 
context of astrophysical grains. 

Measured XAFS spectra for various iron species are shown 
in the Appendix, § 5. The sharp rise in the absorption cross sec- 
tion, also called the edge step, occurs when the incident photon 
is equal to the binding energy Eq of the deep-core electron. The 
cross section, /i(-E), is typically expressed as 

M(£) = Mo(£)[l+x(£)] (1) 

where fJ,o(E) is the cross section of the isolated atom and x(E) 
is the oscillatory modulation of the cross section due to the in- 
teraction of the ejected photoelectron with the condensed state 
material. The features within a few electron volts of the edge 
are called the near-edge structure. The extended X-ray absorp- 
tion fine structure (EXAFS) extends from a few tens to many 
hundreds of eV above the edge. 

The oscillatory portion of Eq. 1, x(E) , is determined by 
the local-scale atomic structure in the vicinity of the absorb- 
ing atom, thus the absorption spectrum is highly sensitive to 
the chemical species of the absorber. This dependence is seen 
in the equation for the oscillatory fine structure in the single 



scattering approximation: 

**> = E ^^ ilm ^R i+ ^m-\ (2) 

In this equation, the summation is over all atoms i surround- 
ing the absorber, Fj(k) and $,(&) are the amplitude and phase of 
the function describing the scattering of a photoelectron from 
neighboring atom ;, R, and of are the distance between absorber 
and scatterer and the thermal RMS variation in that distance, 
is an element-specific, intrinsic loss term (Rehr & Albers 2000), 
X(k) is the mean free path of the photoelectron, and k — the 
wavenumber of the photoelectron — is related to the energy by 

k= y ISf(E-Eo). The mean free path, the l/R 2 dependence, 

and the Debye- Waller-like er 2 term together serve to limit the 
spatial extent probed by the photoelectron. XAFS is, therefore, 
an inherently local probe which measures structure within 10 A 
or less of the absorbing atom. 

As a heuristic explanation of Eq. (2), consider the outwardly 
propagating photoelectron. This propagates as a spherical wave 
and scatters from the surrounding atoms. (See, for example, 
Figure 5 in Rehr & Albers 2000.) This backscattered wave 
interferes with the outwardly propagating photoelectron. The 
density of states of the absorber is modulated by this interfer- 
ence function measured at the position of the absorber. This 
modulation is, according to Eq. (2), a sum of damped sine 
waves in photoelectron wavenumber. 

The frequencies of these sine waves are a measure of the dis- 
tances to the scatterers, as given by the sin(2£7?,) term. The 
Fj(k) in the numerator of the pre-factor depends upon the atomic 
species of the scatterer. The oscillatory fine structure function 
of the EXAFS, is isolated (Newville et al. 1993) by fitting a 
smooth spline representing the atomic contribution to the cross 
section such that the low-frequency Fourier components are re- 
moved from the spectrum. The resulting x(k) can be Fourier 
transformed (Koningsberger & Prins 1988) to produce a com- 
plex function related to the partial pair radial distribution func- 
tion of atoms about the absorber. Examples of the Fourier trans- 
form x(/?) functions are shown for various iron compounds in 
the Appendix, §5. 

With the aid of theory (Zabinsky et al. 1995; Ankudinov et al. 
1998), the x(/?) spectra can be analyzed (Newville 2001) to re- 
veal the details of the local environment about the absorber, 
including the atomic species of and distance to the near neigh- 
bors. Eq. (2) is evaluated for each kind of scatterer in the ma- 
terial. The contributions from the scatterers considered in the 
analysis are summed and then Fourier transformed. The Fourier 
transform of this theoretical x(^) is compared to the Fourier 
transform of the data. The parametric terms in Eq. (2), such as 
R and a 2 , are modified by non-linear minimization techniques 
to best fit the theory to the data (e.g. Newville et al. 1995). The 
Fourier transform is a complex transform, thus this minimiza- 
tion is made by comparing both the real and imaginary parts 
of the transforms of the theoretical and measured spectra. A fit 
to measured data is, therefore, a balance of contributions from 
various scatterers which have subtle amplitude and phase rela- 
tionships. In §5, the magnitudes of the Fourier transforms of 
data from several iron compounds are plotted along with the 
magnitudes of the Fourier transforms of the fitted contributions 
from various near-neighbor atoms. 

Both the near-edge and the extended spectra can be used 
to identify the chemical species of the absorber. The features 
of the near-edge structure of metallic and oxidic iron are dis- 
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tinct. Also the energy of the onset of the oxide edge is shifted 
by ~ lOeV relative to the metal. In the case of the metal- 
lic iron and magnetite comparison considered in §3, the near- 
est neighbor in the metal is another iron atom at a distance of 
about 2.5 A while in magnetite the nearest neighbor is oxygen 
at about 1 .9 A. These differences in distance and species of the 
near neighbor in the various iron species are readily apparent 
in the analysis of the EXAFS spectra, as seen in the Appendix, 
§5. For satellite spectra, we intend to use similar techniques to 
identify the chemical species of grains in the ISM. 

3. CAN WE MEASURE XAFS WITH ASTRO E2 ? 

In order to investigate the limiting effects of noise (§3.2) 
and spectral resolution (§3.3) for detecting an XAFS signal 
from satellite measurements, we have made synchrotron XAFS 
measurements of various Fe compounds, including metallic Fe, 
FeO (wustite), Fe2C>3 (hematite), Fe304 (magnetite), Fe2SiC>4 
(fayalite), and FeS2 (pyrite). In the Appendix (§5; Figs. 3) 
we show the measured absorption spectra of each of these iron 
species along with the spectra convolved by 10 eV to simulate 
the approximate resolution 5 of the Astro E2 instrument. The 
inset to each figure shows the Fourier transform of the isolated 
fine structure spectra. We see that the XAFS data contain am- 
ple information which can be used to distinguish different iron 
species, including the chemical shift of the onset of the edge, 
the features in the near-edge structure, and the details of the 
Fourier transformed fine structure. For the purpose of assessing 
the feasibility of detailed composition studies of dust with As- 
tro E2, we limit our detailed discussion to a comparison of Fe 
metal with magnetite Fe304. Our conclusions, however, are ap- 
plicable to all iron compounds in the vicinity of the Fe K edge 
near 7.1 keV. 

3.1. The synchrotron experiment 

The data presented in this paper were measured at beam- 
line XI 1 A at NSLS or at beamline 13BM at the APS. XI 1A 
is a bending magnet beamline with a double crystal Si(lll) 
monochromator with energy resolution of about 2x1 0~ 4 AE /E. 
Adequate harmonic rejection was accomplished by detuning 
the second crystal to attenuate the incident intensity by about 
30%. 13BM was operated in a substantively similar mode. 
The measured samples were a 12 fim iron foil and a fine pow- 
der of Fe3C>4 dispersed within a quantity of boron nitride and 
cold pressed into a pellet. Both samples were optimized (Lu & 
Stern 1983; Stern & Kim 1981) for high-quality transmission 
mode (Koningsberger & Prins 1988) measurement. The inci- 
dent and transmitted beam intensities were measured by ion- 
ization chambers filled with appropriate mixtures of inert gases 
and operated in a stable voltage regime. 

The synchrotron data were processed using the ATHENA pro- 
gram (Ravel & Newville 2004). Energy calibration was done 
using a reference Fe foil measured in parallel during the mag- 
netite measurement. Normalization was performed in the stan- 
dard fashion (Koningsberger & Prins 1988), i.e. by regressing 
a two-term polynomial to the region before the pre-edge and 
subtracting it from the entire spectrum. A three-term polyno- 
mial was regressed to a region beyond the absorption edge. 
The value of the post-edge polynomial extrapolated back to 
the edge energy was used as the normalization constant. The 
ATHENA program was also used to simulate the effects of noise 
and energy resolution discussed below. Noise was added to the 

5 http://heas arc. gsfc.nasa.gov/docs/astroe/gallery/performance/specres. html 



spectra using a pseudo-random number generator. The RMS 
value of the noise is given as a percentage of the edge step, i.e. 
the normalization constant derived from the polynomial regres- 
sions. Thus a 10% level of noise is understood as having an 
RMS value of 0. 1 on the scale of the normalized spectra shown 
in Figs.l and 2. The effect of resolution was simulated by 
convolving the measured spectra with a Lorentzian of a given 
width. 

3.2. The limiting effects of noise 

At the spectral resolution of the Astro E2 XRS, noise will 
be the most significant impediment to our ability to speciate in- 
terstellar grains. To demonstrate, we investigate the tolerable 
noise level by adding random fluctuations to our synchrotron 
data at the level of 5-10% as described in § 3.1. As seen from 
Fig. 1, a S/N > 10 (or noise <10%) is required to distinguish 
Fe-rich material in gas versus the solid state form to enable a 
direct measurement of the gas-to-dust ratio in a target absorp- 
tion column. At this level, if no XAFS are detected, we would 
conclude that the cloud is predominantly gas phase. 

However, should a given environment be predominantly dust, 
as expected e.g. along the line-of-sight passing through some 
molecular cloud, we would like to be able to distinguish the 
chemical species of the condensed state of the absorber. In the 
comparison of XAFS from iron metal versus an oxidic form 
of iron, Fig. 1 shows that while Fe metal can still be distin- 
guishable at the approximate XRS resolution of 10 eV with a 
relatively high (er rms ^10%) level of noise, a cr rms ~5% noise 
leaves the oxide unrecognizable. Since iron oxides are more 
likely candidates for ISM grains, S/N considerations should be 
one of the most important factors when estimating adequate ex- 
posures times. 

Readers should also take caution, when making analogies 
with astrophysical environments that the numbers given above 
are calculated assuming material in purely solid form. As such, 
when considering astrophysical environments - the dust con- 
tribution to the overall optical depth of the edge should be ac- 
counted for and the numbers given scaled accordingly. To illus- 
trate, the required S/N > 10 value quoted above would translate 
to a required S/N > 33 (i.e. noise <3%) for a pointed obser- 
vation of the diffuse ISM where the dust:to:gas ratio is ~ 30% 
(e.g. Wilms et al. 2000). Molecular clouds should be better tar- 
gets for XAFS studies given the expected higher dust content. 

3.3. The limiting effects of spectral resolution 

We next convolve our synchrotron data of Fe metal and 
Fe3C>4 with Lorentzians to simulate different instrumental reso- 
lutions and find that, to do studies of XAFS at ~ 7.1 keV Fe K, 
a spectral resolution of 10 eV or less (as only currently possible 
with the Astro E2 XRS) is required to extract the X-ray absorp- 
tion fine structure information to deduce the chemical make-up 
of the absorber. As demonstrated in Fig. 2, 20 eV resolution 
with moderate (<7 rms ~5%) noise is still adequate for identify- 
ing metallic iron, but insufficient for distinguishing any kinds 
of oxides. Simulated spectra using the Astro E2 XRS response 
matrix further show that binning the data to increase S/N is un- 
fortunately not an option. A doubling of the S/N would re- 
quire that we bin the data by a factor of 4 (i.e. VAO implying 
a 40 eV resolution - at this resolution, even the high S/N syn- 
chrotron spectra of the different compounds become nearly in- 
distinguishable and vastly more sensitive to noise (Fig. 2 left). 
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4. CONCLUSION 

The unsurpassed spectral resolution of the Astro E2 XRS in 
the > 7.1 keV K-edge region of iron and nickel will allow the 
first detailed study of interstellar grains via analysis of X-ray 
Absorption Fine Structure (XAFS). Given data with sufficiently 
high S/N, the extended portion of the K-edge spectrum can be 
analyzed to reveal information about the chemical species of the 
material in the ISM. As demonstrated here, the spectral resolu- 
tion of Astro E2 will allow us to (1) definitively separate dust 
from gas phase absorption, and (2) give us good expectations 
for being able to directly determine the chemical composition 
of dust grains containing iron. 

We note that, while XAFS theory is quite mature, the study of 
XAFS is still largely empirical. As such, the success of such a 
study will require that space-based measurements (e.g. with the 
Astro-E2 XRS) be compared with empirical XAFS data taken 
at synchrotron beamlines to determine the exact chemical state 



of the astrophysical dust. Because determining the content of 
the ISM has important consequences for understanding the evo- 
lution of the Universe, and X-rays are a powerful tool for en- 
suring the success of such a study, spectral resolution and area 
should remain an important consideration of future X-ray mis- 
sions. 
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FIG. 1 . — The effects of instrument resolution and noise on synchrotron XAFS spectra of Fe metal and Fe304. In both figures, reading from the top, the black 
(solid) lines shows spectra measured at a synchrotron source, the blue (dashed) lines show these data convolved by a Lorentzian of 10 eV width to simulate the 
resolution of the Astro E2 XRS, the subsequent green (dash-dot) line shows the convolved spectrum with artificial noise of ctrms equal to 5% of the step height, and 
the bottom purple (dotted) line shows the convolved spectrum with ctrms equal to 10% of the step height. (LEFT) The unit-normalized absorption cross section /i(£) 
spectra, which emphasize the differences between the metal and the oxide in the near-edge portion of the spectra. The vertical line at 71 12 eV marks the ionization 
energy of metallic iron. For gaseous iron, the oscillatory structure blue-ward of the edge step seen here will be absent. (RIGHT) The magnitude of the complex 
Fourier transform of the X-ray fine structure spectra which identifies the distance to the neighboring atom. The first peak identifies the position and element of the 
neighbor atom. The peak just above 2 A identifies the iron neighbor in metallic iron, whereas the peak below 2 A identifies the oxygen neighbor in an iron oxide. 
Note that convolution and a high level of noise leave the iron metal recognizable as such, while the noise level is a more serious impediment to identifying the oxide 
species which are more likely to be found in the ISM. 




FIG. 2. — (LEFT) The effect of convolution on the Fe and Fe304 spectra. As the width of the Lorentzian convolution is increased, the two spectra become 
increasingly difficult to distinguish. (RIGHT) The effect of noise on the Fourier transformed fine structure becomes more pronounced as the width of the Lorentzian 
convolution increases. At 40 eV convolution width, the peaks near 2 A identifying the chemical species of the absorber become unrecognizable. Therefore, 
increasing the width of the bins is not a suitable method of improving measurement statistics. The full energy resolution of the Astro E2 XRS along with adequate 
time integration is essential. 
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5. APPENDIX: XAFS OF IRON COMPOUNDS 






Energy (eV) Energy (eV) 

FIG. 3. — The effect of convolution on a variety of common iron species: (1) metallic Fe, (2) wiistite FeO, (3) hematite Fe203, (4) magnetite Fe304, (5) fayalite 
Fe2Si04, and (6) pyrite FeS2- The measured synchrotron /i(E) spectra (black lines) are shown near the edge along with those data convolved with a Lorentzian 
of width 10 eV (blue lines) to simulate the approximate energy resolution of the Astro E2 XRS. The vertical line in each figure marks the edge energy of metallic 
iron at 7.112keV. Note that the effects of noise are not considered here. The inset of each picture shows the magnitude of the complex Fourier transform of the 
oscillatory fine structure after it has been isolated from the fi(E) as described in the text for both the raw (black) and the convolved (blue) data. Fits to the raw 
data were performed and the contributions from the near-neighbor iron (cyan), oxygen (red), sulfur (brown), and/or slilicon (green) scatterers are also plotted in the 
insets. Note the variation in scale of the insets — the ability to distinguish by the details of the fine structure in the presence of significant noise is sensitive to the 
spectral strength of the fine structure. 



